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Abstract 

We discuss the ground state and the small-amplitude excitations of dilute boson- 
fermion mixtures confined in spherical harmonic traps at T = 0, assuming repulsive 
boson-boson interactions and with each component being in a single hyperfine state. 
From previous studies of the equilibrium density profiles we propose a schematic 
phase diagram in a plane defined by the variables abf/abb and a^k? , where abb 

and abj are the boson-boson and boson-fermion scattering lengths and k^p is the 
Fermi wave number at the centre of the trap. With this background we turn to the 
equations of motion for density fluctuations in the collisional regime and discuss 
some general features of the eigenmodes. We display analytic solutions for sound 
waves in a quasi-homogeneous mixture and for surface modes at weak fermion-boson 
coupling. 

PACS: 03.75.Fi, 05.30.Jp, 67.40.Db 



1 Introduction 



The achievement of Bose-Einstein condensation in trapped atomic gases [1 - 
3] has given new impulse to the study of dilute quantal fluids. Experimental 
attention has extended to the properties of double Bose condensates [4, 5] and 
more recently to the attainment of quantal degeneracy in gases of fermionic 
atoms [6 - 10] and in dilute mixtures of bosons and fermions [11]. A specific 
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motivation for interest in boson-fermion mixtures is that the collisions between 
fermions and bosons can foster thermalization of the fermionic component 
and induce so-called sympathetic cooling [4, 12]. Elastic collisions between 
fermions in a single hyperfme state are inhibited by the Pauli principle [13] 
and even in multi-component Fermi gases the collisions between atoms in 
different hyperfme states become ineffective at very low temperature because 
of Fermi factors [9] . 

A number of theoretical studies have already been addressed to properties 
of dilute boson-fermion mixtures confined in harmonic traps. They have con- 
cerned the equilibrium density profiles of the two species both at zero and at 
finite temperature [14 - 17], the kinetic energy of the fermionic component [18] 
and collective excitations from sum rules in the collisionless regime [19]. Our 
main focus here is on small-amplitude collective excitations for a mixture un- 
der harmonic confinement in the collisional regime. For a detailed discussion 
on how this regime may be achieved in dilute Fermi gases we refer to the work 
of Amoruso et al. [20]. 

We first set up a general background for our discussion of collective excitations 
in boson-fermion mixtures by constructing a schematic phase diagram from 
the earlier studies of the equilibrium state [14 - 18]. This helps to identify var- 
ious dynamical behaviours, which may later be examined in detail by explicit 
numerical solution of the equations of motion. We then report analytic results 
for the eigenfrequencies of the mixture in limiting situations, corresponding to 
mixing of Bogolubov sound and fermionic first sound in a nearly homogeneous 
mixture and to surface modes in a regime of weak fermion-boson coupling. 

We assume throughout a positive boson-boson scattering length, as is appli- 
cable to Bose condensates of 87 Rb and 23 Na. Mixtures of 6 Li and 7 Li should 
have special interest, once condensation is realized for 7 Li in a hyperfme state 
where the scattering length is positive. According to calculations on the colli- 
sional properties of ultracold potassium [21], the 40 K- 41 K mixture should also 
be an interesting system for study from the present viewpoint, if condensation 
of the rare isotope 41 K can be achieved. 



2 Equilibrium density profiles and phase diagram 

We consider a system of Nf fermions of mass mj and N b bosons of mass m b 
at T = 0, confined by spherically symmetric external potentials Vf ;b (r) = 
m f y bUjj b r 2 1 "2 with frequencies uif and uj b - We assume that a single spin state is 
trapped for each component of the mixture and that all bosons are in the con- 
densate. We can then omit the fermion-fermion interaction, which is inhibited 
by the Pauli exclusion principle. The boson-boson and boson-fermion inter- 
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actions are described by contact potentials with scattering lengths a bb > 
and dbf, yielding coupling-strength parameters g = Aiili 2 a bb /m b > and 
/ = 2%h 2 a b f/m r . where m r is the reduced boson- fermion mass. 

The equations of motion for the partial particle densities n a (r, t) and partial 
current densities j<j(r, t), where a is an index for the atomic species [a — f 
or b) follow from those for the one-body density matrix upon projection on 
the main diagonal (see for instance [22]). They take the form of continuity 
equations, 

d t n a (r,t) = -V-j CT (r,t) (2.1) 

and of generalized hydrodynamic equations which, in a mean-field approxima- 
tion and neglecting coupling to energy fluctuations, can be written as 

m/$j/(r, t) = -V • nW (r, t) - n 7 (r, t) V [V f (r) + fn b (r, t) - H ] (2.2) 

and 

m b d t j b (r,t) = -V-rt b \r,t) 

-n b (r, i)V [V b (r) + gn b (r, t) + fn f (r, t) - fi b ] (2.3) 

In these equations Il^^r, t) are the kinetic stress tensors and \i a are the chem- 
ical potentials of the two atomic species. In the following we treat the kinetic 
stress tensors by a local density approximation for fermions, 

n^(r,0 = ^[n/(r,0] 5 % (2.4) 

with A = h 2 (67c 2 ) 2 / 3 /2rrif and by a Thomas-Fermi approximation for bosons, 

IlW(r,f) = 0. (2.5) 

In these equations we have dropped terms which are quadratic in the velocity 
field. 

The treatment given above assumes that local equilibrium has been established 
for the ground-state density profiles and that it is maintained during dynamic 
fluctuations of the particle densities. The dynamical behaviour that we shall 
discuss in Sect. 3 will therefore correspond to a collisional regime [20]. 
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2. 1 Equilibrium density profiles 



The equilibrium density profiles (n e a q (r), say) follow from Eqs. (2.2) - (2.5) by 
setting the current densities to zero. The results are well known [14, 15]: 

9(pi b -V b (r)-fn e f(r)) (2.6) 
and 

n e /(r) = A~ 3/2 \ N - V f (r) - fn b \r)f 2 0fr f - V f (r) - frtf(r)) . (2.7) 

The chemical potentials are to be determined by self-consistency conditions 
on the particles numbers, N a = J n eq {r)dv. On comparing the two density 
profiles in the limit of vanishing boson-fermion coupling, it is easily seen that 
the quantity (2A/3)[n^ '(r)]^ 1 / 3 may be viewed as an effective fermion-fermion 
coupling arising from the kinetic pressure of the Fermi gas [23]. It is also useful 
to introduce the radii R b = (2/jL b /m b u; b ) 1 ^ 2 and Rf = (2 (if /mfUJ 2 ) 1 ^ 2 for the 
two clouds in the absence of interactions, but taken at the true values of the 
chemical potentials. 

For values of the boson-boson coupling and of the confinement frequencies 
which are typical of current experiments, and assuming comparable mag- 
nitudes for the boson-boson and the boson-fermion coupling, the fermionic 
cloud is considerably more dilute than the bosonic cloud. In this situation the 
boson cloud is essentially unaffected by the interactions with the fermionic 
component and the latter can be treated as confined inside a double-parabola 
effective potential Vf ff (r) [15]. This is given by 

m f uj(l --f)r 2 /2 + f(i b /g forr < R b 

K z -°) 

mjuj 2 r 2 /2 for r > R b 

where 

7 = (fm b u 2 )/(gm f uj 2 ) . (2.9) 

Evidently, if 7 > 1 the effective potential (2.8) has negative concavity inside 
the Bose radius and the fermions will tend to be expelled from the centre of 
the trap. The result at strong coupling is phase separation in the mixture, 
with the fermions forming a shell around the boson cloud [15]. 

A second type of phase separation for the boson-fermion mixture was demon- 
strated in the work of M0lmer [14]. In this case one takes Nf ~ N b and consid- 



n 



cq / 



(r) = g 1 fi b - V b (r) - fn f q {r) 



Vf f {r) = < 
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ers very large values of the boson-boson coupling g, such that the two clouds 
acquire similar densities. One ultimately finds that the bosons are expelled 
from the centre of the trap and form a shell around the fermion cloud. 

Finally, Eq. (2.8) is valid also in the case 7 < 0, i.e. for attractive boson- 
fermion interactions, as long as the fermionic cloud is more dilute than the 
boson cloud. However, the fermions now tend to draw the bosons towards the 
centre of the trap, ultimately leading to collapse once 7 becomes sufficiently 
negative [24]. 



2.2 Schematic phase diagram at zero temperature 



As shown by Vichi et al. [18], the thermodynamics of the boson-fermion mix- 
ture at T = is well characterized by means of two scaling parameters, which 
are 7 in Eq. (2.9) and 



x 




1/2 



V 2m b u b J 
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dbb N b 
a ho (6^)5/6] 



1 1/5 



(2.10) 



with — (h/mbUb) 1 / 2 . The description of the system with only two scal- 
ing parameters, instead of the eight original ones entering the expression of 
the equilibrium density profiles, is a major simplification. The ratio R b /Rf 
determines the deviation of the kinetic energy of the fermion cloud from its 
ideal-gas value [18]. 

For the purposes of the present discussion we make the simplifying assump- 
tions rrif — m b , ujf = uj b and Nf = N b . In this case the two scaling parame- 
ters take the simple expressions x = (15a bb kf ^ / '48) 1 / 5 and 7 = a b f/a bb , with 
kf = {A^N f fl & /a ho being the Fermi wave number at the trap centre for a 
non- interacting Fermi gas. As remarked under Eq. (2.7) above, the quantity 
aff = 37r/4A;j - > may be viewed as an effective fermion-fermion scattering length 
arising from the Pauli kinetic pressure. We can now proceed to construct a 
schematic phase diagram for the boson-fermion mixture in the plane defined 
by the variables 7 = a b f/a bb and y = a bb kf°\ Figure 1 shows the result for the 
case 7 > (i.e. repulsive boson-fermion interactions). 

We search for the boundary of phase separation by means of a simple condition 
of linear stability on the two- by-two matrix of scattering lengths a bb , a b f and 
aff [23] (see also [25]). A mixed state is stable if the inequality a bb aff — a b f > 
holds, i.e. when 
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Fig. 1. Schematic phase diagram of boson- fermion mixtures in a spherical harmonic 
trap, for the case of repulsive boson-boson and boson-fermion interactions. The 
phase-separation region is marked DM (for "demixed"). See the main text for a 
description of areas I, II and III in the miscibility region. 

Within this region we can distinguish a region of appreciable overlap between 
the two clouds (region I in Figure 1), bounded by two regions of diminishing 
overlap: a region at 7 > 1 where the fermions are being expelled from the 
centre of the trap (region II in Figure 1) and by a region at y > 4/3n where 
the bosons are being expelled from the centre of the trap (region III in Figure 
1). If the spherical symmetry of the trapping potential is preserved in the 
demixed state, these two regions end into spherical-shell structures with the 
fermions at large 7 (or the bosons at large y) forming the outer shell. 




3 Small amplitude oscillations around the equilibrium profiles 



The equations of motion for small-amplitude density fluctuations hf(r,t) and 
n b (r, t) in the collisional regime are found by combining Eq. (2.1) with the 
linearized form of Eqs. (2.2) and (2.3). They are 



#n/(r, t) = i-V • [nf{r) V ^A{nJ{r))-^hj{T, t) + fn b (v, t) 



(3.1) 



and 



dtn b (r,t) = — V • W(r)V [gh b (v,t) + /n,(r,f)]} 



(3.2) 
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Again the quantity {2A/3)[nf(r)]- 1 / 3 enters Eq. (3.1) as an effective fermion- 
fermion coupling. 

It is easily checked that the form of Eqs. (3.1) and (3.2) is such as to satisfy 
the generalized Kohn theorem [26, 27] for the centre-of-mass coordinate x(i) 
of the whole fluid when the two confinement frequencies coincide. In this mode 
of motion we have ^&,/(r, t) = n e b q f{r — x(t)) and hence 



for small amplitude oscillations. By substituting Eq. (3.3) in Eqs. (3.1) and 
(3.2) and using Eqs. (2.6) and (2.7), we find 



Therefore, the centre of mass of the system oscillates at the bare trap fre- 
quency. 

Evidently, a general solution of Eqs. (3.1) and (3.2) can only be obtained 
numerically, since the equilibrium density profiles are not known analytically 
from Eqs. (2.6) and (2.7). The dynamics of the fluid as a function of the 
scaling parameters will be especially rich in the regions of the phase diagram 
in which the nature of the partial density profiles is changing on the approach 
to phase separation. Of course, this will be signalled by a softening of the 
modes associated with concentration fluctuations. Below we restrict ourselves 
to analytic results which are applicable to special regimes. 

3.1 Homogeneous limit 

The dynamics of a spatially homogeneous mixture is relevant to the situation 
in which the density profiles of both components are slowly varying in space, 
on the length scale set by the mean interparticle distances. Taking the equilib- 
rium densities as constant, it is easily seen from Eqs. (3.1) and (3.2) that the 
eigenmodes are two sound waves, with linear dispersion relation a>i j2 (<?) = Ci j2 g. 
The sound velocities are given by 



In Eq. (3.5) q, = (gn e ^ /m h ) 1 ' 2 is the speed of Bogolubov sound and Cf = 
[2A(nj (? ) 2 / 3 /3m/] 1 / 2 is that of first sound in a Fermi gas. 



n fci/ (r,t) = -x(t)-V< / (r) 



(3.3) 



m bJ %x(t) • Vn^(r) = JX^V* < / (r)V i V,H,/(r) 



(3.4) 




(3.5) 
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This result may be useful in regard to measurements of the speed of sound in 
elongated systems, in experiments such as already carried out by Matthews et 
al. [30] on a single Bose condensate. 



3.2 Surface modes for weak boson-fermion coupling 



Assuming equal confinement frequencies for the two components of the mix- 
ture (uho, say), it is known from earlier results of Stringari [28] for a Bose 
condensate and of Amoruso et al. [29] for a one-component Fermi gas that 
in the limit / = the frequencies of the surface modes (n= 0) with angular 
momentum number are equal and given by oji = / 1 / 2 cj/ lo . A perturbative treat- 
ment of the effect of fermion-boson coupling is therefore especially simple for 
these modes. 

We define the unperturbed equilibrium density profiles n[°j(r) and the in- 
dependent amplitudes n^'P\r) of the unperturbed density fluctuations, for 
each value of the angular momentum I and of its ^-component m. After 
Fourier-transforming Eqs. (3.1) and (3.2) and linearizing them in the cou- 
pling parameter /, we multiply Eq. (3.1) by [n^(r)] _1 / 3 fij' m ^(r) and Eq. (3.2) 

by n^' m \r), and then integrate both equations over r. We obtain in this way 
a determinantal equation for the shift of the eigenmode frequencies fli due to 
the fermion-boson coupling. Writing 

n? = u 2 ho (l + Ai) , (3.6) 



we find 



a bf 
2a>ho 



(E + B)± siga(a bf )y/(E-B) 2 + 4CD 



(3.7) 



this result being valid for either sign of a b f. The quantities entering Eq. (6) are 
defined in the Appendix, where it is also shown that they can all be expressed 
in terms of the Gauss hypergeometric function [31]. 

The result given in Eq. (3.7) should be better applicable in the case R}y < -R/, 
since in deriving it we have assumed that the equilibrium density profiles for 
both components are perturbed in a symmetric manner. This corresponds to 
working in region I of the phase diagram (and in the corresponding region at 
negative 7). 
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4 Summary and discussion 



The main features of the zero-temperature phase diagram that we have sket- 
ched in Figure 1 for boson-fermion mixtures arise from the competition be- 
tween the kinetic energy of the Fermi gas and the repulsive boson-boson and 
boson-fermion interactions. The former disfavours localization and phase sep- 
aration, while the latter favour the spatial separation of the two components. 
On varying these system parameters while preserving the spherical symme- 
try of the trapped fluid mixture, a spontaneous symmetry breaking occurs 
towards a situation where the two components are spatially separated along 
the radial direction. Evidently, axial separation of the two components would 
instead be observed in an anisotropic trap. One may think of exploiting Fesh- 
bach resonances to tune the scattering lengths towards such phase-separation 
regime. 

Given the complex behaviour of the density profiles on the approach to phase 
separation, a complete investigation of the oscillatory eigenmodes of a boson- 
fermion mixture can only be carried out by numerical means. This will be well 
worth doing once a mixture of specific experimental relevance and with rea- 
sonably known scattering lengths is identified. Of course, the dynamics in the 
phase-separated region is simply related to that of the two pure components, 
aside for the presence of interfacial modes. On the other hand, we have seen 
that suggestive analytic results can be obtained in the quasi-homogeneous 
limit and in the weak boson-coupling regime, for either sign of the boson- 
fermion coupling. With increasingly large attractions between boson and fer- 
mions, however, the mixture is expected to undergo collapse [24]. 
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A Calculation of frequency shifts for surface modes 

We introduce the notations \F) = nj ,m \r) and \B) = n^' m \r), define the op- 
erators P bJ = (//m 6)/ )V-[nJ5(r)V] and Df = (-f /m f )V-{(nf (r)y/ 3 V[(nf (r))- 1 / 3 }}, 

and introduce the scalar-product notation (F\0\F) = J nj (r)Onj (r)dr 
etcetera. With these notations we find the following expressions for the quan- 
tities entering Eq. (3.7): 
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B = (F\(nf\r)r^D f \F)/(F\(nf\r)r^\F) 



Z^r'^m I ixx ™ (x > - x2) ' 1/2(2X ' + n - 3x2UAA) 



C=(F\(nf(r))- 1 ' 3 P f \B)/{F\(nf\r))- 1 /*\F) 

^ r(^/-|-3) r 2/+2 / \r2 2 \ 1/2 { A 0\ 



D={B\P b \F)/{B\B) 

+ (X 2 f - x 2 Y x l 2 \-x 2 + 2(1 + 1)(X 2 - x 2 )}} , (A.3) 

and 



E = -(B\P f \B)/(B\B) 



1(21 + 3) 
2ir 2 X, 



2Z+3 



J dxx 2l+2 (X 2 -x 2 fl 2 . (A.4) 



In these equations T(n) is the Gamma function and we have set X b j = 
Rbj/aho and = dbb/dho- 

All the integrals in Eqs. (A.l) -(A.4) can be expressed through the Gauss 
hypergeometric function 2^1(0, b; c; z), using the relation [31] 
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